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D. A. Garanin and E. M. Chudnovsky 
Physics Department, Lehman College, The City University of New York, 

250 Bedford Park Boulevard West, Bronx, NY 10468-1589, U.S.A. 

(Dated: June 1, 2015) 

Quantum theory of spin relaxation in the elastic environment is revised with account of the 
concept of a phonon spin recently introduced by Zhang and Niu (PRL 2014). Similar to the case of 
the electromagnetic held, the division of the angular momentum associated with elastic deformations 
into the orbital part and the part due to phonon spins proves to be useful for the analysis of the 
balance of the angular momentum. Such analysis sheds important light on microscopic processes 
leading to the Einstein - de Haas effect. 
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I. INTRODUCTION 

The problem of conservation of angular momentum in 
systems containing magnetic moments has been around 
since the discovery of Einstein - de HaaifS and Barnett 2 
effects one hundred years ago. The first effect demon¬ 
strated that the change in the magnetic moment of a feely 
suspended body generates mechanical rotation, while the 
second demonstrated that mechanical rotation induces 
magnetization. For some time the Einstein - de Haas 
and Barnett effects were used to measure the gyromag- 
netic ratio of solids®: The significance of such measure¬ 
ments was diminished by the discovery of the electron 
spin resonance and the ferromagnetic resonance that pro¬ 
vided a more accurate determination of the gyromag- 
netic ratio. After that the experiments on macroscopic 
magneto-mechanical gyroscopic effects have been largely 
abandoned. Surprisingly, however, microscopic mecha¬ 
nisms of the transfer of the spin angular momentum to 
the phonon system and subsequently to the body as a 
whole remain poorly understood. 

The tradition that goes back to the pioneering work on 
spin-phonon relaxation by Van VleclP consists of ignor¬ 
ing conservation of angular momentum under the excuse 
that the Hamiltonian of the system does not possess full 
rotational invariance. It is clear, however, that in theory 
(and in experiment) the angular momentum in a sys¬ 
tem of interacting spins and phonons is conserved. This 
prompted a significant effort by a number of researchers 
to formulate the theory of magneto-elastic interactions 
in a rotationally invariant manner^ 11 ! The advantage 
of such approach is that it is parameter free in a sense 
that spin-phonon rates can be expressed in terms of the 
well-known independently measured parameters. 

Emergence of micro- and nanoelectromechanical de¬ 
vices (MEMS and NEMS) rejuvinated interest to 
the problem of angular momentum in magneto¬ 
mechanical system^®- 2 : Einstein - de Haas effect at the 
nanoscale has been experimentally studied in magnetic 
microcantilever^l^^ and theoretically explained by the 
motion of domain walJ 15 ! Switching of magnetic mo¬ 
ments by mechanical torques in nanocantilevers has been 
proposed^!!!!!. Mechanical resonators containing sin¬ 


gle magnetic molecules have been studied by quantum 
methodifl^H]. Experiments have progressed to the mea¬ 
surement of the angular momentum exchange b etwe en of 
a single molecular spin and a carbon nanotub d 24 * 25 ! 

In nanoresonators the problem is somewhat simpler 
due to the finite number of resonant modes. For a sin¬ 
gle spin in a macroscopic body, however, the number of 
phonon degrees of freedom is practically infinite. In re¬ 
lation to the angular momentum this problem has re¬ 
ceived significant recent at tentio n in experiments with 
atomic spin - based qubits- 2 ® 27 and in application to 
spintroniciPsI To address this problem Zhang and Niu 
recently introduced a concept of the phonon spinPP 

In this paper we investigate this concept for the pro¬ 
cess of the relaxation of a single atomic spin in a macro¬ 
scopic body. By developing an approach similar to that 
for photons we find that within the elastic theory the an¬ 
gular momentum of phonons can be naturally split into 
the orbital angular momentum l/ 1 ) and the spin angular 
momentum lA 2 \ The orbital part corresponds to the ro¬ 
tation of the elastic medium around a certain point, while 
the spin part corresponds to a small-radius circular shear 
displacements of points of the elastic media around their 
equilibrium positions, see Fig. [T] 


Figure 1: Conceptual representation of the motion of the elas¬ 
tic medium that generates the orbital angular momentum L (1) 
and the phonon spin angular momentum L (2) . 

The paper is structured as follows. The concept of the 
angular momentum in classical and quantum theories of 
elasticity is discussed in Section |TT] Conservation of the 












2 


total angular momentum is studied in Section [TTTJby com¬ 
puting its commutator with the Hamiltonian. Quantum 
dynamics of the angular momentum of the relaxing spin 
and emitted phonons is investigated in Section |IV| Sec¬ 
tion |V] contains summary of the results and some final 
comments. 


II. THE ANGULAR MOMENTUM 

A. Angular momentum in the classical theory of 
elasticity 


When a a p is non-symmetric, more care is needed. This 
can happen for two reasons. The first reason is the intrin¬ 
sic anharmonicity of the elastic theory due to the nonlin¬ 
earity of the strain tensor®, 

Uprj = 2 pr] + &r]p T ^-vp^ifrf) (7) 

The fact that H must depend on u pr) leads to 

SH 8u P n SH SH SH 

<j 7<5 — y — 7 7 — 7 f e 7PT (®) 

oe 1 s oe-yg Su pri 8u^$ ou p s 


The angular momentum of the elastic solid is defined 
as 


L = 


d 3 r (r + u) x p, 


(1) 


where time-indendent r corresponds to the non-deformed 
body, u(r, t) is deformation, and p((r, t) = pu((r, t) is the 
momentum density. It consists of two parts 

L = L (1) +L (2) , (2) 

where 

L*' 1 ) = J d 3 rp r x u, = J d 3 rpu x u. (3) 

The orbital part described by L^b corresponds to the 
rotation of the elastic medium around the origin, while 
the spin part described by L^ 2 ) corresponds to a small- 
radius circular shear displacements of points of the elastic 
media around their equilibrium positions, see Fig. [T] 
Applying time derivative to these expressions one ob¬ 
tains 

L (1 ) = J d 3 rpr x ii, f/ 2 ) = J d 3 rpu x ii. (4) 

The dynamical equation for the displacement field is the 
Newton’s equation 


d 2 u a _ da a/3 
^ dt 2 dr@ 


( 5 ) 


with the force in the right-hand-side being a gradient of 
the stress tensor cr a p = SH/5e a p. Here H is the Hamilto¬ 
nian of the system and e a p = du a /drp. After integrating 
by parts in equations Q and assuming zero elastic stress 
at the boundary of the body, one obtains 



d ^a/370'7/33 i-'a') — 




( 6 ) 


In the linear elastic theory i/ 1 ) is zero in the absence 
of internal torques due to the symmetry of the stress 
tensor. Such torques are ignored by the conventional 
elastic theory. L^ 2 ), that is quadratic on deformations, is 
also neglected in the linear elastic theory, making L = 0. 


which is non-symmetric. It is easy to see that in this case 
the second term in Eq. © is needed for the condition 
L = 0 to be exact. 

Anharmonicity, however, is not the only reason for a a p 
to be non-symmetric. It also happens in the presence of 
spins because spin dynamics generates internal torques. 
Consider, e.g., a uniaxial spin Hamiltonian of the form 

H = —D (n • S) 2 (9) 


with n being the magnetic anisotropy axis. Elastic de¬ 
formations of the body rotate the anisotropy axis n by a 
small angle 4> 

0 = iy x u(r), cj) a = ie a/ 3 7 e 7/ 3. (10) 

To the first order in <f> one has n = e z + [(f) x e z ]. Ex¬ 
panding H up to the linear terms in 0, we get H = 
Ha + Hs-ph, where Ha = ~DS 2 and the spin-lattice 
coupling is given bjEl 

-Hs-ph = _ D ( S X S Z + S Z S X ) (j)y + D (Sy S Z + S Z Sy ) (f> X . 

(II) 

The corresponding stress tensor tj a p = SH s - p h/Se a p is 
non-symmetric. Writing it as 

SH s _ ph d-Hg—ph . . 

a< x !3 — c- — 77 7 (lA) 

oe at 3 dcp-r Se a p 

one obtains 

L« = - f ( 13 ) 

J ocf> 

which explicitly expresses the internal mechanical torque 
via the elastic twist. The latter comes from the spin- 
lattice coupling. In what follows we will show that L (2 1 
associated with the phonon spin is also generated in the 
problem of quantum relaxation of the atomic spin. Con¬ 
sequently, L( 2 ), that is usually neglected in the linear 
elastic theory, turns out to be important for the con¬ 
servation of the total angular momentum, even in cases 
when the problem is solved with linear non-interacting 
phonons. 
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B. Quantum theory of phonon angular momentum 

To obtain the second-quantized expression for the 
angular momentum, we use canonical quantization of 
phonons 


, , h e kA e lk r 


h.c 


(14) 


where p is the mass density, V is the volume, e kA are 
polarization vectors, w kA are phonon frequencies and a, 
a' are creation and annihilation operators of phonons. 


One uses Eq. (141, as well as 


p(r) = pu(r) 



e k AA/w k Ae ikr a k A + h.c. 


(15) 

The angular momentum of the body, Eq. (JlJ, consists 
of two contributions, L = LW + i/ 2 ) that have been 
discussed earlier. Here i/ 1 ) is first order in phonon oper¬ 
ators and it can be interpreted as the orbital angular mo¬ 
mentum of the phonons. The term I/ 2 ) is second-order in 
phonon operators and it can be interpreted as the spin 
of the phonons. Splitting the angular momentum into 
two parts is similar to that of photons. It will be shown 
below that the spin of a phonon is h and the phonon-spin 
eigenstates are circularly-polarized phonons. 

The operator of the orbital angular momentum be¬ 
comes 


This operator becomes diagonal in terms of numbers of 
circularly-polarized phonons a k ± = (a k i ± za k 2 ) /a/2 

£ (2) = | (—4 + a k+ + 4_a k _) . (21) 

k 


Each such phonon carries an angular momentum h par¬ 
allel or anti-parallel to its wave vector that can be inter¬ 
preted as the spin of the phonon. 

In what follows we will study conservation of the angu¬ 
lar momentum in the spin-relaxation model by comput¬ 
ing its commutator with the Hamiltonian. Introducing 
spin operators S± = S x ± iS y that follow commutation 
relations [S±,S Z ] = ±S±, one obtains from Eq. ( [TT] ) 

iD 

H a -ph = - 2 is+s z + S Z S +) + h.c., (22) 


where 4>± = (f> x ± i<p y . Using Eqs. (101 and (141 with the 
atomic spin located at r = 0, one obtains 


T>± = 


2 pV 


E 

kA 


e± ■ [ik x e kA ] 
V^kA 


( a kA a kA) 1 (23) 


where e± = e x ± ie y . 


III. CONSERVATION OF THE ANGULAR 
MOMENTUM 

Let us now check conservation of the total angular mo¬ 
mentum 



where j k = i f d 3 r re lkr . As, by symmetry, j k can only 
be directed along k, only transverse phonons contribute 
into i/W In an infinite body, wave vectors are continu¬ 
ous, so that one can replace summation by integration 


J = L + hS 


(24) 


that implies that J must commute with the Hamiltonian. 
The dynamical change of the spin operator has to be 
absorbed by the angular momentum of the elastic matrix, 
whose evolution is given by 


L = 


h L 




ph? 


(25) 


- y\..=> f EE... 

Then one can express j k as 

j k = (27r) 3 9 k <5 (k). 


(17) 


(18) 


Dropping the terms aa and a)a* in L^ 2 ^ that do not 
conserve the number of phonon excitations, one obtains 
after integration over the volume 

L (2) = ^ E t ekA X ekA 'l a kA« k 'A' + h ' C ■ ( 19 ) 

Z kAA' 


In particular, the precession of the spin around the 
anisotropy axis creates the co-wiggling of the elastic ma¬ 
trix with the spin. 

It turns out that by commuting operators one can 
prove conservation of some parts of the angular momen¬ 
tum, whereas the complete prove of conservation requires 
a full quantum-mechanical solution for the relaxing spin 
and phonons created by its precession, presented in the 
next section. The situation is different for the angular 
momentum components perpendicular and parallel to the 
anisotropy axis. We will need commutators 


4>±, L (1) 



E 


(e± • [k x e kA ]) [e kA x j k ] 


(26) 


Keeping only transverse phonons, XX' = 1,2 and using 
[e k i x e k 2] = k /k, one arrives at 

t(2)= ^ET( a k2« k i-4 1 « k 2)- (20) 

k 


<^±,L (2) 



(e± • e kA ) ^a kA - 


( 27 ) 
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that follow from Eqs. (23), (16), and (20). 

Let us first consider dynamics of the transverse com¬ 
ponents of the angular momentum. The dominant source 
of spin precession around the anisotropy axis is the un¬ 
perturbed spin Hamiltonian Ha- 


by a calculation similar to that in Eq. (301. The terms 
* ( 2 ) 

S z and L\ ' are related to spin-lattice relaxation and they 
are sums over k, linear in phonon operators. However, 
one cannot prove 


S x = 


h L 


H a ,S x 


= - l -D [S 2 , 5,] = J ( S Z S V + S y S z ) . 

(28) 

For the matrix, let us first consider the dynamics of the 
phonon orbital angular momentum i/ 1 ). From Eq. (26) 
with the help of the identity 


HS Z + i (2) = 0 


(34) 


Y (e kA -A)(e kA -B) = A-B- (£-A 


A=l,2 


and Eq. (181 one obtains 

h 


= i— Y (e± • [k x e kA ]) (e x 


2V 


B 

(29) 

[e kA x j k ]) 


kA 


Y ( ekA ' t e ± X k ]) ( ekA ' [J k X ^ 


kA 


= I e ± X k l ’ [Jk X e z 


without the full solution of the quantum problem that 
will be presented below. 


IV. QUANTUM THEORY OF THE RELAXING 
SPIN 


A. General solution 

To facilitate solving the problem of spin-lattice relax¬ 
ation, we reduce the spin-phonon Hamiltonian to the 
rotating-wave approximation (RWA) form that conserves 
the energy. Consider transitions of the spin \m — 1) —» 
| to) for to > 0 decreasing its energy and call the spin 
states |1) and 10), respectively. With the help of Eq. (221 
one obtains the spin matrix element of this transition 


= i^y Y {( e± ' jk) ^ e * ' k ) “ ( e ± ' e *) ( k • jk)) 


TO — 1 




s—ph 


%D 

— (2m 1) , (35) 


i\ I d 3 k{k x d k J(k)-(k-d k 6(k))} 


where l m -i m = \JS (S + 1) — to(to — 1). Using Eqs. 
(221 and (|23|, one obtains the RWA coupling in the form 


= ih. 


(30) 


Now from Eqs. (251 and (22) one obtains 


UR 

L„. 


Combining this with Eq. 
tion law 


= -D(SyS Z +S Z Sy). (31) 

, one obtains the conserva- 


where 


kA 


D 




J kA 


^kA — ^ (2TO 1) l-ni—lA 


+ ^kA^ a kA 


H e+ • [k x e kA ] 


2 pV U^kA 


(36) 


(37) 


Hi) 

kiS x T = 0. 


(32) 


r (1) 


L y =0. 


In the same way one can obtain HS y 

However, Eq. (32) is not the whole story. One has 

• (2) in 

to consider L^. „ using Eqs. (271 and (221. The result- 




and the X-operators are defined by 

X ol |l) = |0), X lo |0) = |l). (38) 

The quantum state of the system can be specified by 


ing expression is a sum over k, linear m phonon opera¬ 
tors. It is of the same order as the contribution to hS. 


^ = \cX 


10 


y, c kAO kA ) 100) , 


(39) 


due to the spin-phonon interaction, i 


— ph, Sx t y 


x,y 

that 


kA 


was ignored above. Both terms discussed here are much 
smaller than the dominant terms in the angular momen¬ 
tum, conserved according to Eq. (32 ). These small terms 
are related to the spin-lattice relaxation of the spin. It is 
impossible to prove conservation of these terms without 
performing the full quantum-mechanical solution of the 
problem of spin relaxation. 

Considering dynamics of the longitudinal component 
of the angular momentum, one can prove 


where 100) is the “vacuum” state. T has only one exci¬ 
tation, spin or phonon. Considering the excited state of 
the spin as the reference-energy state, one obtains the 
Schrodinger equation for the coefficients 


5 = -TXykACkA 


kA 


CkA = —i (wkA — W 0 ) c kA - -A kx C, 

n 


(40) 


m = - 

2 h v 


£ s -ph ,4 1} 


= o 


(33) 


where ojq = (Ex — E 0 ) /h is the frequency of the transi¬ 
tion between the spin levels. 






























5 


One can integrate the equations for the phonon modes 
Ck assuming the initial condition of the phonon vacuum: 

i A* r* 

c kX {t) =->c(t') 

h Jo 


a a* r t 

= -/ dre-^^-^cit-T) (41) 

™ J 0 

and insert the result into the equation for the spin c: 

J = ~ 2 E I^a| 2 J* dre-^-^dt - r). (42) 


In this integro-differential equation, c(t — r) is a slow 
function of time, whereas the memory function f{r) = 
SkA |^4kA| 2 is sharply peaked at r = 0. Thus 

one can replace c(t — r) => c(f), after which integration 
over r and keeping only real contribution responsible for 
the relaxation yields the equation 


and thus 


where 


dc r 

dt 2 C ’ 

(43) 

c = e-V^\ 

(44) 

r = ^2 ^kA 2 s (wkA — Wo) 

(45) 


kA 


is the spin relaxation rate. Now, adopting this solution 
in Eq. (411 and integrating over time, one obtains for the 
phonons 


Ck\{t) 


A * e -i(^kA-^oh _ e -(r/ 2 )t 
h w k A - w 0 + *T/2 


(46) 


B. Dynamics of the phonon-spin angular 
momentum 


Let us now compute the phonon-spin angular momen- 
( 2 ) 

turn L\ resulting from the relaxation of the spin. Re¬ 


member that l 1 1} = 0 according to Eq. (|33|) . It is not 


neccessary to use circularly polarized phonons: one can 
work with linearly polarized phonons using Eq. (21 ) and 


Eq. (391. For the quantum expectation value one obtains 


L (2) = ihY^ r ( c k2Cki - 4iCk 2 ) ■ 


(47) 


Using Eq. ( [46] ) and setting WkA => Wk, one obtains 
L (2) = X! (^k2-4ki - 

k 

]_ _|_ e“ rt _ ( e ~ »(“k—wo)t _|_ e i(ui k -ui 0 )t ^ g-(r/2)t 

(w k - w 0 ) 2 + r 2 /4 

( 48 ) 


In the integration over Wk, one goes to the upper and 
lower complex half-plane for the two different oscillating 
terms. As the result one obtains 


L<2) = l 1 _ e rt ) ( Wk ~ w °) ( iA ^ A kl + h.c.). 

(49) 

It remains to show that the integral over k in this expres¬ 
sion can be expressed through T so that T cancels and 
the result simplifies. Indeed, the combination that enters 
Eq. (451 after simplifications becomes 


|2 , i , .2 0 2 iyo— i \ 7 i 2 h . (50) 


|Aki| + |Ak 2 | — D 2 [(2m — 1) 


4 pV O/'k 


On the other hand, in Eq. (491 one obtains 


*A k2 Ak 1 + h.c. = - D 2 [(2m - 1) 


H kk z 
4 pV Wk 


(51) 

Note that in Eq. (|49b only the longitudinal component 

C 2 ) 

L\ is non-zero by symmetry. The latter is just the neg¬ 
ative of Eq. (501 that enters T, Eq. (45). Thus in Eq. 
(49) T cancels out and one obtains the simple behavior 

(52) 


L z = L™ =-(!-■ 


,-rn 


h. 


as the spin undergoes a relaxational transition |m — 1) —» 
|to). This means that the total angular momentum in the 
system spin + phonons is conserved. 


V. DISCUSSION 

We have analyzed the transfer of the angular momen¬ 
tum from the atomic spin to the orbital and spin angular 
momentum of phonons. These two parts of the angular 
momentum of the phonon system are clearly distinguish¬ 
able. The orbital part is first order on the phonon oper¬ 
ators. Its classical counterpart is the twist of the elastic 
matrix around the position of the atomic spin, which is 
linear on the displacement field. The spin part of the 
phonon angular momentum is second order on phonon 
operators. Its classical counterpart corresponds to the 
rotational shear deformations that are quadratic on the 
diplacement field. 

Conservation of the angular momentum in the process 
of the relaxation of the atomic spin has been demon¬ 
strated by us explicitly. It turns out that the change in 
the transverse part of the atomic spin is balanced by the 
orbital part of the phonon angular momentum, while the 
change in the relaxing longitudinal part of the atomic 
spin is balanced by the spin part of the phonon angu¬ 
lar momentum. These findings can be useful in schemes 
where individual atomic spins (e.g., used as qubits) are 
manipulated by phonons. 

An outstanding problem, not addressed in this paper, 
is how the orbital and spin angular momenta carried by 
phonons get transferred to the rotation of the body as a 
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whole in the Einstein - de Haas effect. To answer this 
question one must recall that in a typical Einstein - de 
Haas experiment one induces rotational oscillations of 
a macroscopic body by the low frequency ac magnetic 
field. The corresponding time scales are much greater 
than lifetimes of phonons emitted in atomic spin tran¬ 
sitions. Consequently such phonons fully equilibrate on 
the time scale of the transfer of the angular momentum 


from atomic spins to the body as a whole. 
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